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ABSTRACT 

HCN is becoming a popular choice of molecule for studying star formation in both 
low and high mass regions and for other astrophysical sources from comets to high 
red shift galaxies. However, a major and often overlooked difficulty with HCN is that 
it can exhibit dramatic non-LTE behaviour in its hyperfine line structure. Individual 
hyperfine lines can be strongly boosted or suppressed. In low mass star forming cloud 
observations, this could possibly lead to large errors in the calculation of opacity and 
excitation temperature while in massive star forming clouds, where the hyperfine lines 
are partially blended due to turbulent broadening, errors will arise in infall measure- 
ments that are based on the separation of the peaks in a self-absorbed profile. This is 
because the underlying line shape cannot be known for certain if hyperfine anomalies 
are present. We present a first observational investigation of these anomalies across a 
wide range of conditions and transitions by carrying out a survey of low-mass star- 
less cores (in Taurus and Ophiuchus) and high mass protostellar objects (in the G333 
giant molecular cloud) using hydrogen cyanide (HCN) J=l— >0 and J=3— >2 emission 
lines. We quantify the degree of anomaly in these two rotational levels by considering 
ratios of individual hyperfine lines compared to LTE values. We find that all the cores 
observed demonstrate some degree of anomaly while many of the lines are severely 
anomalous. We conclude that HCN hyperfine anomalies are common in both lines in 
both low mass and high mass protostellar objects and we discuss the differing hy- 
potheses for the generation of the anomalies. In light of the results, we favour a line 
overlap effect for the origins of the anomalies. We discuss the implications for the use 
of HCN as a dynamical tracer and suggest in particular that the J=l— >0, F=0^1 
hyperfine line should be avoided in quantitative calculations. 

Key words: Starless cores - Hyperfine anomalies - ISM: kinematics and dynamics 
- Molecular Clouds - parameter space - dense tracers. 



1 INTRODUCTION 

While hundreds of molecules have now been detected in 
molecular clouds, there are only a handful of molecules that 
are useful as robust dynamical tracers. A good probe species 
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must be abundant enough to be readily observed and it 
must also be excited at high densities, so that it can trace 
the deep dense interior of molecular clouds (where the key 
dynamical processes take place) rather than just the low 
density outer layers. The species must also be chemically 
well behaved so that its abundance relative to hydrogen 
varies predictably. HCN matches all of these qualities and, 
at first sight, appears to be one of the very best tracers 



2 —R. Loughnane et al. 



Low mass J=l-0 



Low mass J=3-2 



M. A 



F=0-l 




-15 -10 -5 5 10 15 

Velocity shift (km s -1 ) 

High mass J=l-0 



12 -10 -8 -6 -4-2 2 

Velocity (km s -1 ) 




-15 -10 -5 5 10 15 



Velocity shift (km s" 1 ) 



A j\ 




70 -65 -60 -55 -50 -45 -40 -35 



Velocity (km s" 1 ) 



Figure 1. Top: Expected line shape for optically thin HCN 
J=l— >0 in a low mass cold quiescent molecular cloud compared 
with a JCMT observation of a low mass protostellar core, L1197. 
The J=l— »0, F=0— »1 hyperfine line, marked 'A' is boosted far 
above its expected strength. Bottom: Expected line shape for op- 
tically thin HCN J=l— »0 in a massive turbulent molecular cloud 
(which increases the line width). This is compared with a Mopra 
observation of a massive core in the G333 cloud. Again the com- 
ponent marked 'A' is boosted but is also much broader than the 
other components. This strongly suggests that line overlap effects 
at higher energy levels are significant. 



of molecular gas in space. HCN was discovered in space 



by Snyder & Buhl ( 1971 1 and examples of the objects it 



has been recently used to observe include comets ||Hoger- 



heijde et al. 2009 Friedel et al. 2005 


Hirota et al.||1999||, 


planetary atmospheres (Marten et al. 


20021, evolved star 


atmospheres (Schilke et al. 2000 Schilke & Menten 20031, 


quiescent nearby low mass stellar nurseries ( Yun et al.|1999 



Park et al.|1999[), dis tant massive star forming regions ( |Hen- 



nemann et al.||2009| |Boonman et al. 20011, active galaxies 



( Kohno et al. |2003 1 and molecular clouds in the high red- 
shift universe (Gao & Solomon 20041. In particular, HCN 



has been used routinely as an infall tracer in low mass star- 
less c ores ( jTafalla et aLpOOo] |Sohn et al.|[2007| |Choi et al 



19991 and also in high mass star forming regions (Wu et al 



2005) 



HCN has a hyperfine structure due to the nuclear 
quadrupole moment of 14 N. This is potentially very useful 
since the optical depth and self- absorption could be mea- 
sured by examining individual hyperfine lines while different 
rotational levels would give a measure of excitation temper- 
ature. 

There is however, a major and often overlooked diffi- 
culty with using HCN for any quantitative calculations in 
that it has a rather pathological and puzzling hyperfine 
structure. In an early paper on this, Walmsley et al. (1982), 
using observations of HCN J=l— >0 in a low mass star form- 
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Figure 2. The hyperfine structure of HCN J=3— »2 compared 
with a JCMT observation of L1622A2. The component marked 
'B' is strongly anomalous. The laboratory spectroscopic rest fre- 
quencies and relative intensities of the hyperfine lines are listed 
in Table Q] 



ing cloud (TMC 1), demonstrated that the hyperfine compo- 
nents are present in ratios that mean they are not in thermal 
equilibrium with each other. The strengths of the individ- 
ual hyperfine lines are 'anomalous' in that they can appear 
boosted or suppressed far beyond what would be expected 
from an LTE analysis (or even from a standard non-LTE 
radiative transfer calculation). This pheno menon of hyper- 



fine anomalies can also be seen in N2H + (Keto & Rybicki 



2010 Daniel et al. 20071, and deuterated nitrogen species 
(Turner 20011 but is not seen in other spe cies with 
hyperfine lines, such as C 17 0, for example (Redman 



et al.||2002j ) 



The underlying mechanism for the HCN anomalies has 
never been fully settled with possible suggested causes in- 
cluding turbulent overlap, radiative scattering and 
line opacity effects (e.g. Guilloteau &: Baudry 1981; 
Cernicharo & Guelin 1987; Gonzalez-Alfonso & Cer- 
nicharo 1993; Turner et al 1997). New observational 
data, coupled with radiative transfer codes that im- 
plement full line overlap offer the possibility of a 
solution to this long standing problem, which will be 
timely in view of the beginning of the Atacama Large Mil- 
limeter Array (ALMA) era. 

In this paper, an attempt is made to systematically in- 
vestigate the HCN anomalies observationally. We present 
the combined results of a HCN molecular line survey car- 
ried out for this purpose toward 25 low mass starless cores 
in Taurus and Ophiuchus and towards 7 massive star form- 
ing turbulent cores in the G333 complex. We find that the 
anomalies are widespread in both classes of object. Our 
study marks the first identification in massive star forming 
regions of the same anomaly as observed in many of the low 
mass cores. We have in addition, for the first time, observed 
two anomalous rotational lines towards the same source. In 
§2 we review the microphysics of the formation of the HCN 
hyperfine lines. §3 details our source selection criteria, the 
molecular transitions and the telescopes used. In §4, the re- 
sults are presented and a method of systematically charac- 
terising the degree of anomalousness is described. An anal- 
ysis of the hyperfine line ratios is then presented. Finally, in 
§5, we summarise our findings and discuss the implications 
of our results. 
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Table 1. Calculated Spectroscopic values for hyperfine compo- 
nents of the J=l— >0, 2— >l and 3— >2 transitions (Aherns 2002) 
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2 THE HCN LINE PROFILE 

HCN has a relatively high dipole moment (/ij = i_).o=2.98E[^] 
for HCN J=1^0 versus ^ J= i^ =0.11D for CO J=1^0) 
so that the lower transitions of HCN, especially, prove to 
be excellent tracers of dense molecular gas in star-forming 
clouds as well as in stellar complexes such as galaxies. This 
is due to the critical densities of rotational transitions obey- 
ing n cr i t oc ^ 2 i/j_ >J _ 1 , for optically thin lines at frequency 
allowing the HCN transition lines to trace ~100-500 



times denser gas than corresponding CO transitions (Pa- 
padopoulos|2007| |. HCN possesses a dominant end nitrogen 



atom that means it is less prone to surface effects on dust 
grain mantles, in particular freeze-out, compared to other 
more abundant H2 mass tracers such as HCO + or C 17 



(Freed & Mangum 2005). Thus HCN remains abundant in 
the gas phase in the cold central regions of star-forming 
cores. 

The (nuclear quadrupole) hyperfine structure of HCN is 
solely due to the non-vanishing electric quadrupole moment 
of the end 14 N nucleus (Qi4 N =20.44 ± 0.03mbQ, which can 
lead to line splitting of the order of several MHz, resulting in 
an easily identified split line structure (with rotational tran- 
sitions in the GHz range). By taking ratios of the relative 
intensities of neighbouring hyperfine lines one can constrain 
the optical depth of a given region in the cloud. Since hyper- 
fine lines belonging to a particular rotational transition differ 
by a few MHz in frequency, there should be little concern 
regarding the coupling of the telescope beam. The physical 
properties of a source such as the density, the temperature, 
and molecular abundance depend on the optical depth and 
so the hyperfine structure should be taken into account upon 
analysing spectral line observations ( Keto fc Rybicki|2010 1. 

Table [l] gives the hyperfine line frequencies and their 
relative intensities. Under LTE or optically thin conditions, 
the ratios of the relative intensities of the neighbouring hy- 



1 1 Debye [D] = 10 -18 statcoul cm 

2 mb => milli-electron-barns, 1 barn=10" 
2 mb milli-clcctron-barns, 1 barn=10" 



perfine lines in each of the two rotational transitions consid- 
ered here generally adhere to a set of fixed constants/limits. 
The relative weightings, given such conditions, in the case 
of the lower rotational transition are of the form 1:5:3 (or 
0.111:0.555:0.333 in terms of relative intensities) for HCN 
3=1— >0. For HCN J=3— >2, four of the six hyperfine lines 
are not spectrally resolved and the spectrum takes the ap- 
pearance of one strong central component and two satellite 
lines, with a ratio of 1:25:1 between them. 

As a preview of the observations presented below, Fig- 
ures [l] and [2] show clear examples of the anomalies. Figure [l] 
shows the theoretically expected optically thin line profile 
for HCN J=l— >0 for turbulent widths appropriate for low 
mass star formation regions (top left panel) and massive star 
formation regions (bottom left panel). Alongside these are 
example observations that clearly show the anomalies: the 
component marked A should go from being the weakest of 
the three hyperfine lines, in the optically thin limit, to equal 
in strength to the other components in the optically thick 
limit. Instead, in both types of region, it is either the broad- 
est or strongest component of the profile. Figure [2] shows 
that a similar effect is present in the higher excitation HCN 
J=3— >2 line. As will be shown in this paper, it is clear that 
the hyperfine anomalies are present not only in the ground 
state of HCN in low mass star forming regions but also in 
other levels and in massive star forming environments. Since 
these anomalies are common this must currently render any 
quantative interpretation of the spectrum of HCN insecure 
at best. 



3 OBSERVATIONS AND RESULTS 

3.1 Low mass star forming clouds 

In order to quantify the presence of HCN anomalies in low 
mass star forming clouds, we collected observations of both 
the HCN J=l->0 and HCN J=3->2 lines for a selection of 
nearby low mass starless cores. 

The sources were primarily chosen from a comprehen- 
sive HCN J=l->0 survey in HCN and HNC, catalogued by 



Sohn et al. ( 2007 I using the Taeduk Radio Astronomical Ob- 
servatory (TRAO) in Korea. The majority of these objects, 
being well-known Taurus- Auriga and Ophiuchus sources, 



were selected by Sohn et al. (2007| as being associated with 
the central peak of the (1,1) inversion transition of NH3, as 
observed by single-dish radio telescopes, and hence tracing 
dense gas ( Lee et al.|1999 1. The objects surveyed are found 
to be at various distances (see Table [6] at the end of the 
manuscript) and have luminosities in the range 0.1 to 0.5 
Lq. ISohn et al. (20071 kindly supplied us with their data 
which form all of our HCN J=l— >0 line profiles. 

The cores to be observed in HCN J=3— >2 were selected 
on the basis of how bright the central F=2— >1 hyperfine com- 
ponent was in the HCN J=l— >0 rotational transition since 
it was anticipated that if this particular hyperfine transi- 
tion was strong then the probability for a detection in the 
J=3— >2 rotational transition of HCN was improved. Cores 
that were also clearly anomalous in the HCN 3=1— >0 line 
were further prioritised. 

Table 2 summarises the telescopes and setups used to 
collect the data. The bulk of the HCN J=3^2 (265.8862 
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Table 2. Observational parameters related to individual telescopes 





Frequency 


Backend 


Au 


A» 




Transition 


(GHz) 




(kHz) 


(kms- 1 ) Fcff 


Vmb 


(1) 


(2) 


(3) 


(4) 


(5) (6) 


(7) 



3- 


=1->Q 


88.6339360 


TRAO auto 


10 


0.033 


0.58 


0.50 


3- 


=3^2 


265.8864343 


JCMT ACSIS 


30.5 


0.034 


0.47 


0.63 


J= 


=3^2 


265.8864200 


JCMT DAS 


78.1 


0.088 


0.47 


0.63 


J= 


=3^2 


265.8864870+ 


IRAM auto 


31 


0.035 


0.68 


0.43 


3- 


=4^3 


354.505500 


KOSMA AOS 


53 


0.059 


0.46 


0.64 



NOTES - (2) fFrequency of the J,F = 3,4— >2, 3 hyperfine component as given by the CDMS, 
http:/ / w ww. ph 1 .uni-ko cln.de / vorhcrsagcn/ 

(6) F e ff — f227r/^4 and f^Tr (H^) is the integral of the beam pattern over the forward hemisphere 
4tt, jBensch et al.|200lf . 



GHz) emission line observations for the low mass sources 
were carried out at the James Clerk Maxwell Telescope 
(JCMT) between September 2007 and July 2008 in band 
5 (t"225 5s0.2) conditions. Single-point observations were ob- 
tained for 30 sources in position-switching mode with a pre- 
defined off-position of 500" in an arbitrary direction. The 
ACSIS digital autocorrelation spectrometer with a band- 
width of 250MHz was used providing a velocity resolution of 
~0. 034km s . The receiver noise temperature (DSB mode) 
was 510-850K. The telescope main-beam efficiency was 0.69 
and the half-power beam width (HPBW) was ~20" in the 
range 225-285GHz. 

Data from supplementary sources were gathered from 
the JCMT archive. This data consisted of HCN J=3->2 line 
observations towards six starless cores L1495A-N, L1544, 
L1517B, L1512, L1521F and L1622A-2 which were carried 
out between February and March 2005, typically in band 
2 (r 2 25 >0.06) conditions. The DAS (which preceded AC- 
SIS at the JCMT) was used in frequency switching mode 
and with a velocity resolution of 0.088 km s . The system 
temperature during the course of this observing run varied 
between 250 and 380K. These starless cores were mapped in 
HCN J=3— >2, in a strip of five successive positions across 
each core. We used the central position data in the analysis 
below but data such as these could be used to investigate 
how the anomalies vary within a single source. 

Additional data, from IRAM and KOSMA, for two 
sources complete our database. These are J=3— >2 and 
J=4— >3 data for the two starless cores, L1498 and TMC-1 
respectively, both part of the Taurus-Auriga complex. The 
former data has been extracted from a molecular survey of 
two prototypical starless cores carried out by Tafalla ct al. 
( 2006 1 whilst the latter core data was supplied from line ob- 

Table|6j at the end 



tation or pulsation that gives doppler shifts in excess of the 
thermal/turbulent line widths (e.g 



servations taken by Ahrens et al 



of the manuscript, lists all the sources observed and gives 
the detections statistics. There were 11 sources with no de- 
tection. 

Figure 3 displays HCN J=l— >0 and J=3— >2 line profiles 
for most of the low mass sources. For the low mass sources 
it can be seen that the HCN J=l— >0 hyperfine lines are well 
separated and resolved. By reference to Figures 1 and 2 it 
is readily apparent that the HCN anomalies are present in 
many of the sources; the leftmost hyperfine line should never 
exceed in strength either of the other two hyperfines. For 
many of the sources, the individual line profiles exhibit ei- 
ther red or blue asymmetry. Such signatures are interpreted 
as being due to dynamical activity such as infall, outflow, ro- 



Mardones et al. 



1997 



Carolan et al.|[2008| [Redman et al.||2004| |Alves, Lada, fc| 



Lada 20011 IRedman, Keto, & Rawlingsl 2006 1 . Examining 



these dynamical properties is beyond the scope of this work 
but many of the low mass sources have been widely investi- 
gated elsewhere. However, it is very much worth pointing out 
that for several of these sources, the HCN J=l— >Q, F=0— >1 
(highest frequency, blue-wards) hyperfine has the reverse of 
the red or blue asymmetry present in the other two hyper- 
fines (e.g. L1521B, L204C-2, L234E-S). This is not likely to 
be explicable in terms of the bulk dynamics taking place in 
the source (see, for example, Stahler & Yen 2010 who model 
such asymmetries in the HCN hyperfine lines, but do not 
address the anomalies) and is instead likely to be part of 
the microphysics of the hyperfine anomalies. This effect is 
returned to in the discussion section. 

As described in §2 the low mass HCN J=3— >2 profile 
can be considered to be composed of a strong central com- 
ponent with two satellite lines. Again, the anomalies are 
readily apparent in many of the sources. The left satellite 
line is heavily suppressed and the right satellite is boosted 
to comparable or higher strength than the central compo- 
nent. The central component is marginally resolved at the 
thermal and turbulent line widths present in low mass star 
forming regions so the line shape is complex as a result of 
partial blending and dynamical effects in the source. The 
right satellite would provide a better guide to the underly- 
ing line shape of each hyperfine in these conditions. For the 
sources with detections, Table [7] at the end of manuscript, 
gives the antenna temperature of the strongest hyperfine and 
the integrated HCN J=l-K) and HCN J=3->2 hyperfine line 
intensities. 



3.2 Massive star forming regions 

An excellent example of a massive star forming region is 
the G333 molecular cloud associated with the RCW 106 
Hn region. This giant molecular cloud complex spans an 
~ 0.7 deg 2 region, and is ~ 3.6 kpc away. G333 has been the 
subject of a large multi-molecular line legacy survey using 
the 22-m Mopra Telescope of the Australia Telescope Na- 
tional F acility (ATNF) between 2004 and September 200 6 
(see e.g. ( |Lo et al.|200"9||Wong et al.|2008| [Bains et al.|2006[ )). 
The data presented in this paper were collected with the 8- 
GHz wide band Mopra spectrometer (MOPS) centred at 87 
GHz. The FWHM beamsize is «36" ( |Ladd et al.|2005| ). Full 
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Table 3. HCN J=l— >0 Observations towards G333 cores 



Emission 




R.A. 


Dec 






v lsr 


Peak 


(J2000) 


(J2000) 


(K) 


(kms- 1 ) 


(1) 




(2) 


(3) 




(4) 


(5) 


G332. 741-0.619 (MMS73) 


16 


20 07.0 


-50 59 


50 


0.864 


-50.28 


G332.772-0.589 (MMS71) 


16 


20 05.7 


-50 57 


04 


1.372 


-55.94 


G332.813-0.703 (MMS75) 


16 


20 49.7 


-51 00 


15 


0.755 


-53.61 


G333.096-0.503 (MMS50) 


16 


21 14.2 


-50 39 


11 


1.861 


-55.57 


G333.241-0.517 (MMS38) 


16 


21 52.5 


-50 34 


13 


0.740 


-49.09 


G333.293-0.423 (MMS30) 


16 


21 42.9 


-50 28 


09 


1.440 


-50.18 


G333.297-0.357 (MMS26) 


16 


21 25.6 


-50 24 


15 


1.527 


-49.82 



NOTES.- (1) - Name of source using its galactic coordinates and Mookerjea 
et al (2004) identifier; (2). (3) - Right Ascension (R.A.) and Declination 
(Dec.) of each source in their corresponding equatorial coordinates; (4) An- 
tenna temperature of the central F=2— »1 hyperfine component; (5) Local 
standard of rest (LSR) velocity of central hyperfine component. 



details of the observational setup is available in |Lo et ah] 
(20091. 



We used this Mopra data set to identify individual mas- 
sive star forming cores in HCN J=l— >0 emission. Because 
this region is distant and very confused, and multiplicity is 
common in massive star formation, it was important for the 
present study to identify clear isolated examples of massive 
cores. HNC is often co-spatial to HCN but is less abun- 
dant (though see Sarrasin et al 2010) and exhibits no hy- 
perfine structur^] Therefore, clumps of HNC emission that 
exhibit single peaked line profiles are likely signposts of iso- 
lated dense cores that will be bright in HCN. We identified 
seven new HCN sources in the G333 data cube in this way. 
These sources were detected in the 1.2 mm emission cata- 
logue G333 of Mookerjea et al (2004) and their cold tem- 
peratures and dust masses ranging from ~ 400 — 1200 Mq, 
mark them as cold massive cores. The positions of the seven 
sources are displayed in Figure 4 and are noted in Table 3. 

Figure [5] displays the HCN J=l— >0 line profiles from 
these seven cores. For these massive core HCN J=l— >0 lines, 
the three individual hyperfine components are still clearly 
identifiable but, due to the greater degree of turbulence, 
the line widths are much greater and are comparable to the 
separation of the hyperfines. The resulting hyperfine line 
parameters for the HCN J=l— >0 data for each of the seven 
massive cores are summarised in Table 3. The anomalies 
here are of a notably different character in that now the line 
widths of the individual hyperfines can be different as well 
as the line strengths. The presence of a variation in the line 
widths of individual hyperfine lines within a single rotational 
transition is an important result in the context of the origin 
of the anomalies and is returned to in the discussion section 
later in the paper. 



4 HYPERFINE INTENSITY ANALYSIS 

Reviewing Figures 1 and 2 and the full catalogue of line pro- 
files in Figures 3 and 5, the anomalies are readily apparent 



3 HNC does actually possess an underlying hyperfine structure, 
but the location of the 14 N-nucleus, means that the individual 
components are unresolvable and a single peak is observed 



visually. In order to characterise the anomalies in a quan- 
tative fashion, we here define the methods used to quantify 
the degree of anomalousness in the two rotational transi- 
tions. The character of the anomalies is then analysed for 
the two types of star forming region. 



4.1 HCN J=1^0 

For the HCN J=l— >0 transitions, we use the method first 



utilised by Gottlieb et al. ( 1975 1 where the ratios of the in- 



dividual hyperfine component integrated intensities, I, are 
measured and compared with Local Thermodynamic Equi- 
librium (LTE) values. These ratios, R(I[F=1->1]/I[F=2->1]) 
and R(I[F=0->1]/I[F=2->1]) (hereafter, simply R 12 and 
R02) were calculated for each core and compared with the 
ranges expected for non-anomalous values. 

The ratios R02 and R12 equal 0.2 and 0.6, respec- 
tively, under LTE or optically thin conditions. Under op- 
tically thick conditions, the lines tend to saturation such 
that they are of equal intensity i.e. R02 = R01 = 1. Thus 
the expected, non- anomalous, ranges are 0.2<Ro2<l and 
0.6<Ri2<l. Cores with hyperfine spectra displaying ratios 
lying outside each of the two ranges are described as be- 
ing the most anomalous, with those outside only one of 
the above ranges being intermediary anomalous candidates. 
Those sources of emission with ratios lying within the two 
ranges represent the non-anomalous cores. The integrated 
HCN J=l— >0 line intensity ratios are recorded in Table 4. 



4.2 HCN J=3^2 

There are six possible hyperfine lines in HCN J=3— >2 but 
due to blending (the hyperfine splitting decreases strongly 
with rotational level) it is not possible to resolve the cen- 
tral group of hyperfine lines and the spectrum takes the ap- 
pearance of one strong central component and two satellite 
lines. The central group of three hyperfines can be identi- 
fied as a collective branch that result from a net change of 
1 in the total angular momentum quantum number, AF=1. 
One hyperfine branch also belonging to this central group 
is negligibly weak. The remaining two branches each result 
in an isolated hyperfine line. With this blending, the op- 
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Figure 4. HCO + J=l— >0 (left) and 13 CO J=l— >0 (right) integrated intensity (zeroth-moment) maps of the G333 giant molecular cloud. 
The crosses mark the positions of the seven HCN anomaly emission sources presented in this work. 
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Figure 5. HCN J=l— >0 line profiles from the seven G333 GMC sources. 



tically thin form of the intensity weightings is 1:25:1 (or 
0.037:[0.200+0.296+0.428]:0.037 =>■ 0.037:0.925:0.037). 

It is practical then to measure the ratio of the outlying 
hyperfines to the central component in order to quantify 
the degree of anomalousness. We introduce the expressions 
R 



AF(0- 



.!), hereafter R -i and Ro+i which 
represent the ratios of the central hyperfine branch (AF = 1) 
compared with that of leftmost (AF = 0~) and rightmost 
(AF = + ) hyperfine branches respectively. The integrated 
HCN J=3— >2 line intensity ratios are recorded in Table 4. 



4.3 Anomalies in low mass star forming regions 

Table 4 summarises our analysis of the collected datasets of 
the two transitions towards the twenty-nine low mass star- 
less cores that were detected in both lines. Examination of 
the table reveals that while the anomalies are common, the 
degree of anomalousness varies between sources. For exam- 
ple, the lower transition ratios are Ro2=0.95 and Ri2=0.80 
for L1197, and R 02 = 1.18 and Ri 2 =1.08 for L1544 (recall 
that a value of 1 in the case of either ratio represents the 
LTE condition of saturation in each of the respective emis- 
sion lines). In the higher transition, the values of the ratios 
are R -i = 0.28 (six times in excess of the LTE value) and 
R 0+1 = 0.78 in the case of L1197 while for L1544 the ratios 
are R - 1 — 0.64 (fifteen times in excess of the expected LTE 
value) and Ro+i = 0.58. Thus it can be seen that L1544 is 
more anomalous than its Ophiuchus counterpart. 



Figure 6 displays N 2 H + J=l— >0 column density calcu- 
lated from optical depth data from Crapsi et aL| |2005), who 
observed the majority of our sources, against the R02 (top) 
and R12 (bottom) ratios. Recall that in LTE these both 
of these line ratios must be in the range 0.2 < R02 < 1 
and 0.6 < R12 < 1. Each of the hyperfine intensity ratios, 
R02 and R12 demonstrate a steady increase as a function 
of the column density of an optically thin tracing species 
(N 2 H + J = 1 — ¥ 0). This result is largely in keeping with 



earlier predictions (Zinchcnko & Pirogov 1987). The cor- 



relation coefficients calculated from the plotted points are 
r = 0.5282 for R 02 and r = 0.5216 for R 12 . We consider 
these to be strong correlation coefficients and therefore con- 
clude that the denser the optically thin line along the line of 
sight, the more anomalous the corresponding HCN J=l— >0 
spectrum and that, from this data, it appears that it is the 
R12 hyperfine ratio that is most responsive to the N 2 H + 
J=l— >0 column density. 



4.4 Anomalies in massive star forming regions 

As noted above, a crucial finding in our G333 massive cores 
is that in all cases one or both of the two outlying HCN 
J=l— >0 hyperfines are significantly broader than the cen- 
tral hyperfine. For example, the line widths of the three 
hyperfine components of G333. 297-0. 357 are 7.8, 4.2 and 
3.4 km s _1 , a variation of a factor of more than two. As 
a result of this, to quantify the anomalies, we use the in- 
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Table 4. Hypcrfinc intensity ratios in the HCN J=l— >0 
and J=3— »2 transitions assembled for each of the 29 
detected cores in the higher transition. 
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a The ratios R„_^ and R.Q+ ^ are defined in relation to the 
J — 3— >2 hypcrfinc structural branches in §4.2. 



tegrated line intensity rather than just the line strength. 
This method, similar to that used by |Lapinov ( |1989[ ) is 
described in the Appendix. The hyperfine intensity ratios 
towards high mass cores were formally computed on the 
basis of the ratio of their respective line strengths, e.g. 
Roi=Ta(F = -> 1)/Ta(F = 2 1). This approach has 
bee n app lied to the IR source NGC 7538 IRS1 by |Cao et| 
al. (19931, where the respective hyperfine components were 



subject to considerable partial convergence on account of dif- 
ferent broadening schemes and microturbulence. The results 
of the analysis of the G333 HCN J=l— »0 data are presented 
in Table 5 which demonstrates that while the majority of 
these cores are anomalous in this transition, the line ratios 
do not exceed unity. 

We do not have corresponding HCN J=3— >2 data for 
the G333 sources but line profiles have been obtained in this 
transition towards some massive star forming regions. Ex- 
amples can be seen in Wu & Evans (2003), Wu et al (2005) 
and Carolan et al (2009). At first glance, it would be ex- 
pected that the hyperfine structure is not resolved in these 
sources. However, Carolan et al (2009) found that in order 
to correctly model the HCN J=3— >2 spectra, the underly- 
ing hyperfine line structure needed to be taken into account 
in order to correctly measure dynamical processes such as 
infall. Regarding any hyperfine anomalies, the heavy blend- 



ing would render an analysis similar to that undertaken for 
the HCN J=3— >2 low mass source data very difficult if not 
impossible. 



4.5 Characterising the anomalies 

In Figure [7] we plot the R02 ratio against the R12 ratio for 
both sets of cores. For the low-mass cores, this involved plot- 
ting the ratios after calculations using the integrated inten- 
sity values in Table [jj with the respective errors along each 
dimension also computed. In relation to the higher mass 
cores, the G333 objects were plotted using the ratios calcu- 
lated from the values indicated Table 5. In the high-mass 
plot, Fig. [7|b), we also include the corresponding hyperfine 
intensity ratios for the eight sample HCN J=l— >0 spectra for 



the sources indicated in figure 1 of Pirogov ( 1999 1. For both 



sets of readings, we computed corresponding errors based 
on the analysis using the technique outlined in Appendix A. 
It can be seen from comparing Figs.[7|a) and[7|b) that the 
ratios in the high mass cores are more tightly constrained 
with respect to their collective hyperfine ratios than the low 
mass sources. This, together with the fact that the high 
mass cores considered as part of this study do not display 
hyperfine ratios exceeding unity, indicate that warmer re- 
gions tend to form anomalous HCN J=l— >0 spectra that 
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Table 5. Summary of HCN J=l— >0 Hyperfine Component 
Analysis towards G333 Cores 
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NOTES. - (1) - Source name in galactic coordinates; (2) - 
Respective hyperfine component; (3) - Line width of each 
hyperfine component in kms" 1 ; (4) - Integrated intensity 
of individual hyperfine components in K km s — 1 . 



favour the pumping of the central F=2— >1 hyperfine com- 
ponent, that is, the thermal overlap description considered 
in Guilloteau & Baudry (19811. This explanation will be 



returned to in §5. Due to the relative low densities of our 
candidate low-mass objects, these figures may also reveal 
that the conditions conducive to such regions favour path- 
ways of excitation at the quantum level that differ from the 
LTE routes, as discussed in 55. 



5 DISCUSSION 

5.1 Hypotheses for the origin of HCN hyperfine 
anomalies 

The HCN hyperfine anomalies are unlikely to be due to 
a simple non-LTE opacity effect because i) the anomalies 
are seen in regions of varying column density and non- 
anomalous clouds are seen with a similar range of den- 
sity ( Lapinov||1989 1 and ii) non-LTE radiative transfer cal- 
culations do not reproduce the anomalies ( [Stabler fc Yen| 
2010). Several hypotheses have been presented in the last 



few decades to account for the anomalies. Radiative trapping 
in the hyperfine transitions ( |Kwan &; Scoville|1975| ), can be 
ruled out upon closer examination of the molecular physics 
involved ( LapinovflT9 89). Collisional rates for HCN for 
conditions in dark clouds are now well established 
(see Sarrasin et al. 2010 for the latest rates, not- 
ing also Faure et al. 2007) so poor rate data are 
unlikely to be the cause of the anomalies. Gottlieb et 
al. ( 1975 1 developed an idea which attempted to address the 



problem by way of probable de-excited photons from the hy- 
perfine lines of upper rotational transitions. This introduced 
the possibility of line overlap for the first time. Guilloteau| 
|fc Bau dry (1981) developed on this theory in their thermal 
overlap treatment of the 3=2— >1 transition for clouds where 
T=30K. Their treatment demonstrated an overpopulation 
of the J=1,F=2 level at the expense of the J=1,F=1 and 
F=0 levels, respectively. With regard to the J=l— >0 rota- 
tional transition, such an emphasis in population leads to a 
strengthening of the F=2— >1 main hyperfine line and a cor- 
responding weakening the two side hyperfine components. 
However, our and other observations towards low mass cores 
show a relative strengthening in these side components so a 
simple thermal overlap effect involving a single level can be 
discounted. 



5.2 A line overlap effect for the HCN hyperfine 
anomalies 

The above studies did demonstrate that it is possible to 
change the strength of individual hyperfine lines through 
line overlap and this was further investigated by Zinchencko 
& Pirogov (1987). A line overlap treatment involving sev- 
eral rotational levels and including line broadening due to a 
combination of thermal, doppler and turbulent effects could 
then offer a promising possible explanation for the anoma- 
lies. There are several pieces of evidence from our study 
that supports such an explanation. The lowest HCN transi- 
tion gives rise to lines that are widely spaced in frequency. 
Transitions between higher rotational levels lead to reduced 
separations between the hyperfine lines, leading to blending 
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Figure 6. Hyperfine ratios R02 (top) and R12 (bottom) versus N2H+1 
the top left corner. 



>0 column densities. The correlation coefficient r is recorded in 



beyond around the HCN J=4— >3 level in low mass clouds 
and the HCN J=3— >2 level in high mass clouds. The hy- 
perfine lines arising from transitions between the J=4,3 and 
2 levels could then be sensitive to doppler effects arising 
from, for example, enhanced or suppressed stimulated emis- 
sion due to locally doppler shifted photons. This is perhaps 
what is been seen in two effects noted earlier in the paper. 
Firstly, in our low mass star forming clouds, several sources 
exhibited a mixture of blue and red asymmetric line pro- 
files in the same rotational level transition but in all three 
cases (and noting the very small number statistics) it was 
the leftmost line (J=l— >0, F=0— >1) that differed from the 
other two. Secondly, in our high mass star forming clouds, 
the hyperfine line widths varied within the same rotational 
level transition. It is the leftmost hyperfine line involved 
again, this time being typically significantly wider than the 
other lines. A line overlap in a higher rotational level tran- 
sition could be the cause of these effects. In the very narrow 
lines in the low mass cloud, a small change in line strength 
at certain wavelengths could disturb the line shape. In the 
broad lines in the high mass cloud, the leftmost transition 



appears to have grown in width at the expense of the central 
component in particular. 

In the line overlap hypothesis, the anomalies in- 
crease with optical depth, which may also explain 
the reason why an isotopomer of HCN, H 13 CN, 
does not exhibit anomalous intensity ratios (e.g. for 
TMC-1 [Irvine & Schloerb|1984| ). This species, which 
is chemically identical to HCN and with an identi- 
cal dipole moment, is lower in abundance by a factor 
of 40-50 and so is therefore typically optically thin. 
The apparent absence of an anomalous feature in 
this tracers hyperfine spectrum and its close relation 
to the more abundant HCN, indicates a role opac- 
ity effects must play in the formation of anomalous 
intensities. 

A limitation of our observational approach is 
that for the majority of sources, we were constrained 
to a single pointing towards the centre of the source. 
Cernicharo et al (1984) presented the first large- 
scale map of the HCN J=l— >0 anomalies, in TMC-1. 
The anomalies were seen to be spatially extended 
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Figure 7. HCN J=l— >0 hyperfine ratios, R12 plotted against R02, for (a) the 64 low-mass starless sources in Table 1 of |Sohn et alTj 
2007] ! and (b) the 1 5 high -mass objects considered in this work, made up from the seven G333 cores indicated in Fig. |??| and the eight 
sources in figure 1 of Pirogov (1999). In (c), both sets of object's ratios are plotted in tandem. LTE ratios for different optical thicknesses 
would lie on the dashed curves in each figure. 



but at the low opacity envelope at the edge of 
the source, the central F=2^T component remains 
strong relative to the two side components (|Guil-| 
loteau &j Baudry||1981[ ICernicharo & Guelin||1987l) 



while closer to the centre, the F=0^T component 
is at its strongest. This could be another manifes- 
tation of the opacity dependence for the anomalies 
discussed above. Cernicharo et al (1984) describe a 
concept by which radiative scattering of the photons 
emitted from the high density core is taking place 
at the lower density periphery. This mechanism can 
reproduce the extremely small HCO + /H 13 CO + inten- 



sity ratios observed towards Heiles Cloud 2 (which 
contains TMC-1) ( |Cernicharo 8c Gue"Iin||1987| ) and 
could play a significant role in the radiative transfer 
of HCN. 

Another possible contributing factor to the presence of 
the anomalies involves considering the allowed downward 
pathways of excitation to the lower rotational state hyper- 
fine energy levels. We examined all possible downward path- 
ways from the J— 7 level to the J=0 level in HCN. There are 
36 different allowed routes via the different hyperfine tran- 
sitions from this upper level, J=7, to the lower J=l level. 
Remarkably, there is only one pathway out of the 36 that 
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leads to emission of the leftmost lower hyperne transition in 
HCN, the J=l— >0, F=0— >1 line (which is involved most in 
the two effects described above). Three of the J=3-2 com- 
ponent hyperne lines also only have one allowed route to 
them. When we considered the same situation for N2H + , 
the number of pathways leading to line excitations as a re- 
sult of transitions to each of the lower hyperfine levels, JFiF 
= Oil and 012, is 2-2.5 times that leading to excitation due 
to transitions to the JFiF = 010 level, thereby exhibiting a 
similar pattern to HCN above. In contrast, a similar ex- 
ercise for a non-anomalous hyperfine species, C 17 0, 
shows an even number of routes down to the lower state hy- 
perne energy levels. Therefore, it seems reasonable to sup- 
pose that small changes to the transition probabilities due 
to line overlap will readily disturb excitation of the HCN 
J=l— ?>0, F=0— >1 line in particular as well as to transitions 
down to the JFiF=010 hyperfine level in the case of N2H + . 
The C 17 line profile is much more robust because the lines 
are assembled from multiple downward hyperfine excitations 
in an homologous manner, that is, the hyperfine lines in the 
lower rotational state are excited equally. This may be due 
to a steady and consistent decrease of the available hyperfine 
energy levels as the hyperfine energy ladder is dismounted. 
Thus, if the number of pathways of excitation coming from 
upper hyperfine energy states and going to lower hyperfine 
energy states from a given hyperfine state (defined as the ref- 
erence hyperfine energy state) does not diminish regularly 
as a result of a particular framework of hyperfine energy 
states, this will lead to a magnification of preferred routes 
of excitation down to the J=l hyperfine levels (as is seen in 
HCN and N2H + ), perturbations of which result in hyperfine 
anomalies. It is important to note that this is a simplified 
argument that only addresses the downward excitation of 
available hyperfine energy states governed by the selection 
rules. 

It is worth noting that N2H + is also reported to exhibit 
anomalous hyperfine spectral components towards some 
sources ( Keto fc Rybicki|2010 Daniel et al.|20 07) but not in 
G333. Since N2H" 1 " is found to be concentrate in the centre of 
dense quiescent cores and not participating in the large scale 
dynamical phenomenon traced in G333 by species such as 
HCN, this again points to a velocity doppler shifting causing 
line overlap being the likely cause of the anomalies. 

5.3 Observational implications 

Our study has shown that the hyperfine anomalies are com- 
mon in star forming regions. If the line overlap hypothesis 
discussed above is correct, what are the implications for the 
use of HCN data in such regions? 

Clearly, measuring opacities by examining the relative 
intensities of the HCN hyperfine lines alone is unsafe and 
alternative methods should be used. HCN can still be used 
cautiously as a dynamical tracer species however. In regions 
where the turbulent velocity is low enough to resolve the 
HCN J=l— >-0 hyperfine structure and discern self- absorption 
effects, then since the F=0— »1 appears to be most easily dis- 
turbed by line overlap effects, it will present a line profile 
shape that could be a complicated mixture of dynamical and 
selective absorption or emission effects. The line width may 
be affected too. Thus this line is an unreliable dynamical 
tracer in any source exhibiting resolved HCN (1-0) hyper- 



fine anomalies and the central component should be used 
instead. In more turbulent sources such as in massive star 
formation, the J=l— HI, F=0— >1 hyperfine will give an in- 
accurate measurement of the line width (but of course a 
lower abundance species without hyperfine structure would 
be selected for such a measurement). In summary, it would 
seem prudent to simply exclude the HCN J=l— »0, F=0— >1 
hyperfine line from any quantative calculation. 

For HCN J=3— >2 for low mass, low turbulent width 
sources where the hyperfine lines are partially resolved, the 
central component should be reliable as a dynamical tracer 
since it is composed of three blended lines with many possi- 
ble downward transitions to them. However, since the left- 
most hyperfine is typically suppressed and the rightmost 
satellite boosted, then in circumstances where the lines are 
not quite resolved the velocity centroid may be miscalcu- 
lated as lying between the central peak and rightmost satel- 
lite giving an error of order ~ 1 km s _1 (half the separation 
of the central peak and rightmost satellite). In massive star 
forming regions, the possible presence of the anomalies mean 
that it should not be assumed that this line is composed of 
the central component alone. As discussed in section 4.4, the 
underlying hyperfine structure must be included to attempt 
to model infall and outflow using HCN J=3->2 in massive 
star forming regions (Carolan et al 2009) and if the anoma- 
lies are present in the way found for low mass sources then 
the line shape will be distorted. In particular, constraining 
the degree of infall by measurements of the degree of line 
splitting (such as via the method of Mardones et al 1997; 
see Wu & Evans 2003 for an example) is likely to result in 
an overestimate compared with using a species like HCO + . 
The magnitude of the difference in infall measurements be- 
tween HCN and an alternative tracer species will be of order 
~ 2 km s" 1 , the separation of the central peak and right- 
most satellite. 



5.4 Radiative transfer and future work 

A full detailed radiative transfer calculation of the line over- 
lap hypothesis is clearly required to verify if the qualita- 
tive explanation above is correct. Line overlap, quite of- 
ten termed line fluorescence has been successfully used at 
other wavelengths to explain, for example He II optical line 



emission from ionized nebulae ( Kastner & Bhatia 1990 \ 



and X-ray Fe lines very close to accreting massive black 
holes ( |Fabian et al.| [2000). Line overlap is very difficult 
to calculate numerically and is typically not avail- 
able in publically released radiative transfer codes. 
However, individual studies have succeeded in us- 
ing line overlap to model the hyperfine line anoma- 
lies. Gonzalez-Alfonso & Cernicharo (1993) used a 
modified Monte Carlo method to analyse the hyper- 
fine anomalies in terms of the scattering approach 
described above. Daniel et al. (2007) have analysed 
the hyperfine component intensities of N2H + for sev- 
eral clouds by deriving a consistent set of hyperfine 
collisional coefficients for this molecule ( |Daniel et 
[aI71|2005D as well as carrying out nonlocal radiative 
transfer calculations of observational data using a 



large velocity gradient model (Daniel et al. 2006). 
Recently, |Daniel &: Cernicharo| ( |2008[ ) implemented 
the Gauss-Seidel algorithm in spherical geometry 
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and included the case of line overlap of hyperfine 
transitions. 

As part of our future work, we are modifying 
our RT co de MOLLIE to use the ( jRybicki &: Hum-| 



1992) line overlap algorithm in order to model 



the HCN anomalies. Our intention is to carry out 
a constraining of the parameter space giving rise to 
hyperfine anomalies in each of the two rotational 
transitions tabulated in Table [7] It is hoped that the 
physical conditions common to our full contingent of twenty- 
eight starless cores will allow for such a study to be done at 
two individual transitions, while our database of 65 starless 
cores in the HCN J=l— >0 transition will allow for the study 
to be carried out extensively at the lower transition. Such a 
project will attempt to identify the physical characteristics 
of the star-forming core that gives rise to the anomalous in- 
tensities. This will focus on physical parameters such as the 
density, thermal as well as turbulent contributions to the 
line width, the velocity field of the core as well as the tem- 
perature gradient. Such constraining will involve the mod- 
elling of each core in multiple transitions of several different 
species with a 3D molecular radiative transport code ( |Keto| 
et al.|2004||Redman et"aI1|2002[ ). 
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6 CONCLUSIONS 

While the existence of HCN anomalies have long been recog- 



nised in the J=l— 5>0 line (e.g. for TMC1, Walmsley et al. 
1982 ) for low mass star forming sources, our survey shows 



for the first time that the J=3— >2 can also be dramatically 
anomalous (see Figs. 1 and 2) and that these anomalies are 
common. Our study shows that massive star forming re- 
gions also exhibit hyperfine anomalies in the J=l— >0 line. 
These anomalies, in contrast to those of the low mass star 
forming sources, are apparent in line widths rather than line 
strengths. 

We favour a line overlap effect for the origin of the 
anomalies. It is likely that in higher rotational levels, where 
the hyperfine lines are more closely bunched, line overlap 
leads to preferred and suppressed radiative decay routes 
down to the lowest energy levels, where the hyperfines are 
widely separated, therefore, emphasising the disproportion- 
ate intensities. An attraction of the line overlap hypothesis 
is that is not strongly dependent on cloud column density 
since it is the form of the thermal, velocity and turbulent 
widths in the cloud that are the trigger for the anomalies. 
In a future paper we will carry out a full radiative transfer 
calculation of the hyperfine spectrum of HCN. 
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APPENDIX A: CONVERGED HCN J=l->0 
HYPERFINE COMPONENT ANALYSIS 

The technique that was employed to deal with the partially 
overlapped hyperfine components for the J=l— >0 transition 
is outlined as follows. Initially, for each component in an 
overlapped hyperfine spectrum, a gaussian approximation 
was implemented (using Starlink Spectral Analysis Tool, 
splat) that took the following form: 



y(x) = Aexp -0.5 



(Al) 



where A is the scale height (peak height) of the approxi- 



mated region (the gaussian peaks in Fig. Al I, c is the central 
velocity-position of the gaussian-peak and a is the gaussian 
width. Thus once the fitted gaussian function is satisfactory, 
there will be a unique set of values for A, c and a, based on 
the form of the original spectral component. 

In our proposed exact calculation of the overlapped re- 
gion, the area is equivalent to the sum of two contributions. 
Firstly, in Fig. |A1| the section of the common area to the 
right of the line marked by v m id (v m id marks the velocity- 
position where the intensity of the overlapped area reaches 
a maximum) is equal to the integrated area over the com- 
plete gaussian on the LHS minus the integral over this gaus- 
sian function, yi(v), from the limits -co to v m id. Secondly, it 
is reasonably straightforward to deduce that the remaining 
section of the "common area" to the left of the line marked 
by v m id is equal to the integrated area over the complete 
gaussian on the RHS minus the integral over this gaussian 
function, y2(v), from the limits v m id to +00. The overall 
area, E, is then the sum of these two segments giving, 



E = 



+ 



yi(v)dv - / yi(v)dv 

) J —OO 

r- 00 

y2(v)dv — / y 2 (v)dv 

> J v mid 



(A2) 



The integration over a gaussian function above or be- 
low a certain threshold v m id is not available analytically, but 
the solution has a simple relationship to the error function, 
Erf(v), or its complement, Erfc(v) (tabulated in mathemat- 
ical tables, Beyer 1991), 



yi(v) = ya(v) 

After much rearranging of Eq. 
quadratic in v is achieved, 



(A4) 

the following 



/ 2 

(0-2 
1a\a\ In 



0"x)v 2 + 2((TiC2 



2 \ 
C2C1) 



2 2 
+ °lC 2 



2 2 

C"2Cl 







(A5) 



Using the standard techniques for solving Eq. ( A5 1, two 
values are achieved, one of which lies between the values ci 
and C2 each for the two gaussians respectively. This is v m id in 
Eq. (A2| above. The other value is another velocity-position 
at which the two overlapping gaussians intercept, but far 
from the vicinity of "common area" in Fig. |A1| given above 
(this depends on the width of each of the respective gaussian 
functions). Neglecting the latter, the equation for E above 
can be solved routinely to give the area of the overlapped 
region. This area is apportioned to each of the overlapped 
spectral components (these spectral components are actu- 
ally hyperfine branches) with respect to the ratio of the non- 
overlapping sections of these profiles, i.e. for the LHS com- 
ponent in Fig. |A1| above, the non-overlapping area is given 
by f!™" d yi(v)dv. The middle hyperfine branch in the HCN 
J=3— >2 transition overlaps both the left and right branches, 
and so the non-overlapping segment of this branch was de- 
termined by subtracting sections from each of the overlapped 
areas. Once the respective portions from each of the over- 
lapped regions were distributed accordingly to each of the 
hyperfine components, the relative intensity ratios, R02 and 
R12 could be calculated. Fig. [7] includes these ratios, each 
value subject to the analytical treatment described here. 



Erf(v) 
Erfc(v) 



u du 



V 71 " ■ 
1 - Erf(v) 
2 



du 



(A3a) 



(A3b) 



It is the error function itself that plays an important 
role in each of the above integrals in Eq. \A2\ . Additionally, 
Erf(0)=0 and Erf(oo) = l aided in the calculations. Along 
with the unique set of constants that describe the form of 
the fitted gaussian function for each of the hyperfine spectral 
components, the quantity, v m id, crucial for the integrated 



emission calculation via Eq. ( A2 1, had to be determined an- 



alytically for each set of overlapping gaussian components. 
In accordance with the above description, each of the gaus- 
sian approximated functions are represented by yi(v) and 
ya(v), respectively. Then, via Eq. \A1\ , the calculation pro- 
ceeded as follows: 
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Figure Al. Resultant gaussian fit to the J=l— >0 hypcrfine components of the bright infrared source 06056+2131 (source data from 
|Pirogov|1999| ). 



